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NOISE REDUCTION IN IMAGES 
FIELD OF INVENTION 

The present invention relates generally to noise reduction, 
and specifically to noise reduction in digital images. 

5 BACKGROUND 

All known methods of reduction of high frequency, random noise 
in images reduce to replacement of initial pixel value by 
some average value over some area near the initial pixel. 
They differ mainly in the choice of filters, that is in the 
10 choice of the weights, with which the pixel values are summed 
during averaging^- and in the choice of the summation 
algorithms saving computing time. The summation algorithms 
may involve ^ Fourier transformations, but the principle 
remains the same. 

15 These methods do not distinguish the noise from the useful 
high-frequency components of the true image and as a 
consequence they reduce both to the same extent. The noise 
reduction smears an image, and the fine details of the image, 
if they were obscured by the noise, become even less visible 

20 after the noise reduction. The smearing is ' an unavoidable 
consequence of noise reduction, if no additional information 
about the noise or the image structure is available. 

SUMMARY OF THE INVENTION 

An object of the present invention is to improve noise 
25 reduction in images, and to diminish negative consequences of 
noise reduction . 

This object is among others attained by methods as claimed in 
claim 1, 6, or 9, and by a computer program product as 
claimed in claim 16. 
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An advantage of the present invention is achieved by 
utilization of different methods in dependence on information 
known of noise or image structure, whereby improved noise 
reduction is achieved in accordance with said information. 

5 Further features and advantages of the present invention will 
be evident from the following description. 

DETAILED DESCRIPTION OF EMBODIMENTS 

In the following description, for purpose of explanation and 
not limitation, specific details are set forth, such as 

10 particular techniques and applications in order to provide a 
thorough understanding of the present invention. However, it 
will be apparent -to one skilled in the art that the present 
invention may be practiced in other embodiment that depart 
from these specific details. In other instances, detailed 

15 descriptions of well-known methods and apparatuses are 
omitted so as not to obscure the description of the present 
invention with unnecessary details. 

The present invention can, if some additional information 
about noise or image structure, besides the image itself, is 
20 given, be used to improve noise reduction and reduce the 
image-degrading consequence incurred by the noise reduction. 
Four important examples of such additional information is the 
knowledge that an image has one or several of the following 
properties : 

25 (a) almost all the area of the true image, i.e. of the noise 
free image, is covered by many-pixel fragments, wherein the 
intensity is changing smoothly; 

(b) the dispersion or other parameters of noise reduction, 
such as FWHM (full width at half maximum) , are approximately 
30 known; 
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(c) the noise in the image is significantly correlated with 
the noise in other images (or in different parts of the same 
image) ; and 

(d) the image is computed (reconstructed) from source data, 
5 which are common to source data of other images. 

Most images have property (a) , though the sizes and shapes of 
smooth fragments may vary much from image to image. 

Property (b) is valid for most noisy images, including X-ray 
pictures . 

10 Properties (c) and (d) are valid for most reconstructed 
images, including tomographic images and maps of bone density 
reconstructed from dual-energy X-ray images. 

Although the noise in an image may be correlated with the 
noise in other images or portions thereof for many technical 

15 reasons, the most frequent cause of such correlation is the 
origin of images from the same source data. When several 
images are computed (reconstructed) from the same data, the 
same initial noise enters into all the computed images 
thereby making the noise, obtained by reconstruction, 

20 correlated. 

The noise in reconstructed images is usually strongly 
amplified compared to the noise in source data and becomes an 
imminent problem, which can be solved by the methods of the 
present invention as will be, described below. An important 

25 example is the calculation of bone and soft tissue densities 
from two X-ray images of the same part of body, but made with 
X-rays of different energies. Since the reconstruction of the 
densities is based on the small differences between such X- 
ray images, the reconstructed densities are very sensitive to 

30 the noise in these source images and the noises in the 
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reconstructed images are relatively large and are strongly 
correlated with each other. 

The properties (a), (b) , (c) , and (d) give new possibilities 
of noise reduction which were not fully exploited before. 

5 Methods of noise reduction will be described which using 
properties (a), (b) , (c) , and (d) accomplish deeper noise 
reduction and diminish the smearing consequence incurred by 
it. Especially deep noise reduction without any smearing of 
small details is achievable for the reconstructed images when 

10 methods using properties (a), (b) , (c) , and (d) are combined. 

Noise in an image is generally any contamination dirt 
covering the true -image of the object of interest. The shapes 
of dirt fragments and the statistical properties of dirt may 
be quite different. The present invention concerns mainly the 

15 noise, which looks as dark and light spots or grains 
scattered over the whole image and covering it densely. If 
the image of the same or different object is produced again, 
e.g. another X-ray image taken of a patient, these spots lie 
at different places in that image than in the first image. 

20 The high-frequency noise in pixel images consists of small 
spots of the size of one or two pixels. Such noise is often 
seen in images made by high sensitive films and electronic 
cameras made in conditions of poor illumination or low X-ray 
flux, when the number of registered photons over a pixel area 

25 is not large. 

The pixel value p at some point X in a noisy image can be 
considered as a sum of a mean value P, which the pixel would 
have in a true image, without noise, and of fluctuation F, 
wherein the pixel value p is given by the formula: 



30 



p = P 4- F. 
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The expression "mean value" is denoted by the symbol M, so 
that P = M(p), and M(F) = 0. The main parameter of such 
noise, describing its strength, is dispersion D. The 
dispersion D(X) is defined as a mean value of the squared 
5 fluctuation for pixel X, which is given by the formula: 

D(X) = M(F 2 ) . 

The quality of an image is usually characterized by the so- 
called signal-to-noise ratio (S/N) , which is actually defined 
as P 2 /D. The noise reduction makes the dispersion D smaller 
10 and improves S/N. 

The mean value M(FG) of the product of fluctuations F, G of 
two pixel values -p - P(X) + F, g = P(Y) + G at different 
points X, Y (Y may. be the point with same coordinates as X 
but lying in a different image) is called covariance of 
15 values p, g. This quantity together with dispersions of p, g 
makes 2x2 covariance matrix 

\^21 ^22 J 

with elements, 

V n = M(F 2 ) - D(X) , 

20 V12 = V 2i = M(FG) , 

V 22 = M (G 2 ) - D(Y) . 

The strength of the mutual dependence of fluctuations F, G is 
measured by the correlation coefficient C, which is given by 
the formula: 



25 



C = V 12 /VVn^ 
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which may vary from -1 to 1. The values C = -1 and C - 1 
correspond to complete dependence of fluctuations, so F and G 
are proportional and by knowing F one may thereby calculate 
G, and vice versa. The dispersion D of the linear combination 
<xp + Pg of the pixel values p, g is expressed by the formula: 

D(ap + pg) = a 2 Vu + p 2 V 22 + 2apV 12 . 

The mean value is most often unknown and is approximately 
estimated by the average value over some number N of 
independently obtained values x± of the same quantity x, 
wherein the average value of x is given by the formula: 

Avr (x) - X^^XX ' 

where positive coefficients w± are weights. If the weights 
are equal to one, which is adopted if x ± have same accuracy, 
the average of x is given by the formula: 

Avr(x) - ^x ± /N . 

For example, the true pixel value P(X) at point X may be 
estimated as Avr(p) of values p± of pixels at points Y± in 
some area around X, if there are reasons to think that true 
values Pi differ little from P(X), or their deviations from 
P(X) compensate each other. The dispersion D(X) may be 
estimated as the average of the squared deviations from 
Avr(p), which is given by the formula: 

D(x) « Avr((p - Avr(p)) 2 ) , 

if there are reasons to think that the noise is of the same 
strength in the chosen area around X. Similarly, covariance 
V12 may be estimated as: 

V 12 * Avr((p-Avr(p)Xg-Avr(g))) . 
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The usual methods of noise reduction replace each pixel of an 
image by the average value over some area around a pixel. 
This area together with a table of coefficients, with which 
the pixel values are summed, is usually called a filter. A 

5 simple example of a filter is a square of 9 pixels, the 
values of which are summed with the same coefficient 1/9. when 
the average is calculated. Averaging with such filter 
diminishes the noise dispersion 9 times, but smears the image 
making each visible line in it 4 pixels wider and 

10 correspondingly decreasing its contrast. The methods 
described in the present invention achieve greater reduction 
of noise dispersion at the same smearing consequence, or 
reduce the smearing consequence for the same noise reduction. 

A tomographic reconstruction is the calculation of the 
15 density of the object on some planes cutting the object using 
a plurality of X-ray pictures or NMR (Nuclear Magnetic 
Resonance) data. It is a complex procedure, reconstructing 
images with correlated noise. 

The reconstruction of the densities of the bone and of the 
20 soft tissue is more simple and consists of calculation of two 
functions a(p,g), b(p,g) giving the densities of the bone and 
soft tissue as functions of pixel values p,g at point X in 
two X-ray images obtained with X-rays of different energies. 
The functions a(p,g), b(p,g) are usually obtained from 
25 calibration X-ray exposures of objects with known densities 
of bone and soft tissue equivalents. Both functions are very 
sensitive to small differences of pixel values p - g, so the 
calculated densities have larger noise dispersions and worse 
signal-to-noise ratio than the original X-ray images. The 
30 noise in the bone density image is strongly correlated with 
the noise in the soft tissue density image, where the 
correlation coefficient C is negative and close to -1. 
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A first method of noise reduction using property (a) of the 
true image (presence of smooth fragments) and property (b) of 
the noise (known dispersion) will now be described. 

The first method reduces noise at each point X of an image by 
averaging the pixel values over a region R, wherein the 
averaging region R is dependent on X and is selected as the 
largest region around X, which includes only pairs of pixels 
Y, Z symmetrically placed with respect to X and having the 
half sums (p(Y) + p(Z))/2 of values deviating from each other 
within the limits corresponding to the level of the noise in 
the image and to a user-defined tolerance preferences L for 
the noise-reduced image. 

The largest averaging region R approximating the optimal one 
is found as follows'. 

Firstly, a selected pixel X is included into the region R. 
Then, the pairs of pixels Y=X-v, Z = X + v, where v is a 
shift vector, symmetrical with respect to point X and 
touching the already filled part of the region R, are 
considered pair by pair and a check is made whether the 
squared difference of their half sum p(v) = (p(Y) + p(Z))/2 
from the average value of the pixels already included into R 
does not exceed the dispersion D of the noise present in the 
mentioned difference multiplied by some tolerance level L set 
by the user. If the pair of pixels passes the test, it is 
included into R. As long as the region R growths, the process 
continues with considering of" pairs as described above. The 
process is stopped when no new pair of pixels passes the 
tests. Then the average pixel value over R is used as a value 
of pixel X in an image with reduced noise. 

Better approximations of an optimal region R may be found in 
a similar way by inclusion of more tests for randomness of 
deviations of pixel values and tests for the presence near X 



of certain expected non-random details of the true image like 
thin lines. The found region R may be revised and made more 
dense and regular by inclusion into the region R of part of 
the internal points rejected in the first pass and by 
smoothing the borderline of the region R. 

The region R must be symmetric to eliminate the distortion of 
the image due to contribution of the gradient of true image 
intensity into the average value. 

The limitation of deviations of half sums of pixel values 
from the average value over R limits the distortion of the 
noise-reduced image due to the contribution of the curvature 
of the true image ^intensity by about Vdl/4 . 

Generally, the region R' is elongated in the direction of 
minimal curvature of the image intensity. For pixels X near 
the boundary of smooth fragments, the region R becomes narrow 
and may be reduced to a line. 

Due to statistical uncertainty of the boundary between the 
region R and surrounding image fragments, the averaging over 
the region R does some smearing of the image. However, the 
range of such smearing is several times smaller than the 
range of the smearing in the usual noise filtering methods, 
which is of the order of the filter diameter. 

The choice of the averaging region R symmetric with respect 
to X is usual for all filters making the noise reduction, 
since it makes the average closer to the true value than in 
case of an asymmetric region. The present invention discloses 
that the region R is optimized individually for each point X 
and that in the selection of this region only the half sums 
p(v) of pixel values at symmetrical points Y, Z are used 
instead of individual pixel values. This makes the process of 
choice of region R insensitive to the gradient of intensity 
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of the true image. In particular, this property has the 
above-mentioned test written as the formula: 

(p(v) - Avr v (X) ) 2 < D(X,v)L, 

where Avr v here denotes the average over all pixel pairs 
5 earlier included into R and D(X,v) stands for the dispersion 
of the difference p(v) - Avr v (X). 

Dispersion D(X,v) - D(X)*(l/2 + 1/n) , where n is a number of 
pixels in R, may be known from different sources: from 
analysis of earlier images obtained with the same technology, 

10 from the number of photons registered in case of Poisson 
statistics of noise, or from the image itself by means of the 
formula D{X,v) ='"Avr((p - Avr(p)) 2 ), where p = p (v) and 
averaging is done oyer all pixel pairs earlier included in to 
R. The estimation of D(X) and D(X,v) may be done in many ways 

15 and is not an integral part of the invention. 

The details of the algorithm picking up the pairs of pixels 
tested for possible inclusion into region R may be different 
and are not the integral part of the present invention. The 
process of checking is mentioned solely to give a general 

20 idea how the method may be realized. Actually, in a tested 
algorithm, the order, in which the pairs of pixels were 
tested, was fixed so as to allow the region R to grow 
continuously, by one-pixel layers of square shape until the 
deviation checks were positive. Every time when a new pair 

25 was added to the region R, the pixels touching to this pair 
were marked as possible extensions of R. The next picked pair 
of pixels was first checked for being marked. If it was 
unmarked, it was skipped without other more time-consuming 
tests. When none of the pixels in a next completed square 

30 layer passed the tests, the process of expansion of the 
region R stops. 
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The purpose of the additional checks for randomness of 
deviations is intended to clear the region R from compact 
groups of pixels (and their symmetrical partners) which pass 
the tests for individual pairs, but which deviate 
collectively too much to one side. If deviations are random, 
the squared deviation of the average of K half sums of pixel 
pairs from the true value (P(X) - Avr K p(v)) 2 may only exceed 
2D(X)/K in a few percent of cases. If a group of pixel pairs 
deviates more, its exclusion from the region R will most 
probably make the average over the region R closer to the 
pixel value P(X) of the true image. The check for the 
presence of certain suspected details is similar to the check 
above, but is more sensitive for deviation in groups of 
pixels of certain shapes. 

The region R found in a one-pass process usually contains 
many holes, where most of the pixels do not deviate much, but 
at the time of their check during the first pass did not 
touch pixels already included into R. The revision of R with 
the account of points included afterwards and with slightly 
more liberal deviation test does not increase the size of R 
but makes it more solid. 

The pixels in a borderline of a solid many-pixel region R 
should, normally, touch with 4 or 5 pixels belonging to this 
region. The rejection of peripheral points touching with less 
than 3 pixels belonging to the region R does not decrease 
significantly the number of pixels in the region R, but 
smoothes the borderline and decreases the risk of smearing 
distortions . 

A second method of noise reduction using property (c) that 
the noise in the image is significantly correlated with the 
noise in other images will now be described. 
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The second method reduces noise at each point X of an image 
Ii correlated with the noise at a corresponding point Y of an 
image I 2 with the variance-covariance matrix V for pixel 
values p(X,Ii), p(Y,I 2 ) by: 

5 (1) obtaining from image I 2 a noise reduced image I 3 by any 
noise reduction method not using image Ii; 

(2) calculating the estimate F(Y) given by the formula: 

F(Y) - p(Y,I 2 ) - p(Y,I 3 ) 
of the fluctuation of the pixel value at point Y; 
10 (3) calculating the estimate H(X,Ii) given by the formula: 

„ ( x,i 1 , . Ew y'^^'ji 

v(p(Y,I 2 ),p(Y,I 2 )) 

of fluctuation G(X,Ii) of pixel value p(X,Ii) in image Ii; and 

(4) obtaining from image Ii a noise reduced image l 4 by 
subtracting from each pixel value of image Ii the estimate of 
15 its fluctuation, which is given by the formula: 

p(X,I 4 ) = p(X,l!) - H(X f Ii) . 

If the matrix V as a function of X is known, implementation 
is straightforward and reduces to repeating calculations (2), 

(3) , and (4) for all pixels X. 

20 If the matrix V is not given,- its estimate at each point X 
may be found by a standard method as average values over a 
small region around the point X and corresponding region 
around the point Y of products of deviations of pixel values 
from their average values over the same regions. 

25 The second method reduces in image I x only the part of its 
noise which is correlated with the noise in image I 2 and does 
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not remove the part of noise in image I 1 which is independent 
from the noise in image I 2 . So the maximal coefficient of 
noise reduction by the second method is limited by the 
difference 1 - C 2 , where C is the correlation coefficient, 
which is given by the formula: 

c ^ Vfe^ljpfo,!,)) 

A /v(p(x,I 1 >p(x / l 1 ))v(p(Y / I 2 ) / p(Y / I 2 )) * 

If C is close to 1.0, the noise reduction in image I 4 may be 
very deep. If C 2 is noticeably smaller than 1.0, image I 4 
still contains noticeable independent noise. This noise may 
be reduced by applying to image I 4 other noise-reduction 
methods (including the first method) . 

The noise reduction by the second method has some specific 
consequence. It is related with the fact that image I 3 
obtained from image I 2 contains smearing distortions. These 
distortions pass into estimates H and from them pass into the 
noise-reduced image I 4 , where they look as a kind of shadows 
of boundaries of objects on image I 2 . 

Since images I lF I 2 are different, the boundaries of image 
fragments in image I 1 and image I 2 , generally do not coincide 
and shadows of boundaries of fragments of image I 2 do not 
smear details in image I 4 . These shadows obscure details in 
image I 4 much less than the original noise in image Ii, so 
small details in image I x lost in the noise in image I x 
become visible after such noise reduction. By contrast, usual 
noise-reduction methods make small details less visible after 
noise reduction . 

The intensity and widths of mentioned shadows of image I 2 in 
image I 4 depend on the choice of the method and parameters of 
noise reduction for image I 3 . The optimal choice of noise 
reduction for image I 2 depends on the structure of images Ii 
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and I 2 and of preferences of the user of image I 4 . The choice 
of the first method of the present invention for noise 
reduction in image I 2 is advantageous since it reduces the 
shadows of image I 2 obscuring the image I 4 . 

The noise remaining in image I 4 can be further reduced by 
many methods. A comparison of image, obtained from image I 4 
after noise reduction by the first method, with a true clean 
image shows that, when the noise correlation is present, the 
combination of the second and the first method makes very 
deep noise reduction without noticeable distortions or 
smearing . 

The roles of images Ii and I 2 may be interchanged, so the 
noise may be reduced in image I 2 by the second method as 
well. 

In an experiment with two artificially constructed images Ii, 
I 2 , imitating images taken with high energy respective low 
energy X-rays, the noise was 95% correlated. To reduce noise 
in image I 2 , the first method was used with an additional 
check sensitive to the systematic deviation of local 
intensity from the average one and reducing "shadows" of 
image I 2 in image I 4 . 

A comparison of image I 4 and noise-reduced image, obtained 
from image Ii by the first method, gave the result that the 
first method is more efficient for large fragments (for 
background) , while the second method is advantageous for 
small fragments and details. 

A third method of noise reduction using properties (a) , (b) , 
(c) , and property (d) that images are reconstructed from 
common source images (data) will now be described. 

The third method reduces noise in images I x and I 2 
reconstructed (computed) as known smooth functions Ii = 
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a (81,82), I2 = b(Si,S 2 ) from source images Si, S 2 with 
independent noise of known dispersion by combination of the 
first and the second method modified to exploit images Si, S 2 
and functions a, b as follows: 

(5) a noise reduced image I 3 is obtained from image I 2 by 
noise reduction by the first method, where the averaging 
region R is selected as the maximal region R around Y, which 
includes only the pairs of pixels symmetrically placed with 
respect to Y and; either 

having the sums of values deviating from each other within 
the limits corresponding to the level of the noise in the 
image and to a tolerance preferences of the user of the 
noise-reduced image; or 

corresponding to pairs of pixels in source images Si, S 2 
having the sums of values deviating from each other within 
the limits corresponding to the level of the noise in the 
images Si, S 2 and to a tolerance preferences of the user of 
the noise-reduced image; 

(6) the estimate F(Y) given by the formula: 

F(Y) = p(Y,I 2 ) - P(Y,I 3 ) 

of the fluctuation of the pixel value at point Y is 
calculated; 

(7) the covariance matrix V for pixel values p(X,Ii), p(Y,I 2 ) 
at points X, Y on images Ii, I 2 is computed in linear 
approximation as given by the formulae: 

V(p(x,l 1 ),p(x,l 1 )) = | J^-l D, +(#-1 D 2/ 







D 1 + 




) 


l 5S 2 J 
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v(p(x,i 1 ), p (y,i 2 ))= 



f a_ V a,_ A 



da 



db 



fda Y*^ 



D 1 + 



V 9S 2 j\v^ 2 j 



db 

as, 



D 2 , 



v( p (y,i 2 ), p (y,i 2 ))= 



'Sb^ 2 



v 5s iy 



v 5S 2y 



D„ 



where Di and D 2 are the noise dispersions in images Si, S2 at 
points Z, T corresponding to point Y; 

5 (8) the estimate H(X,Ii) given by the formula: 



H(X,Ii) = f(y) 



vCpCx^J^pCy,!,)) 
v( p (y,i 2 ) /P (y,i 2 )) 



of fluctuation G(X, Ii) of pixel value p(X,Ii) in image Ii is 
calculated; and 

(9) a noise reduced image I4 is obtained from image Ii by 
10 subtracting the estimate of fluctuation from each pixel 
value, wherein P(X,I 4 ) is given by the formula: 

p(X,I 4 ) = p(X,Ii) - H(X,Ii) . 

The relations Ii = a(Si,S 2 ) and I 2 = b(Si,S 2 ) between images 
are pixel-wise, that is the point X on image Ii corresponds 
15 uniquely to points Y, Z, T on images I 2 , Si, S 2 respectively. 
So, these relations can be written for pixels as given by the 
formulae : 

p(X,Ix) = a(p(Z,Si) ,p(T,S 2 ) ) , 

p(Y,I 2 ) = b(p(Z,S!) ,p(T,S 2 ) ) . 

20 The implementation is straightforward and consists of 
performing operations (5) -(9) for all pixels of image Ii. The 
simplest implementation of point (5) is to check only the 
difference between the sum of values and the double' value of 
a central pixel (the check is done both for image Ii and for 
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the source images Si, S 2 ) ♦ The use of source images in the 
step (5) usually makes the quality of the region R found by 
the pixel-selection process described in the simplest 
implementation of the first method good enough without 
additional checks and boundary smoothing. 

More advanced implementation of point (5) may check as well 
the deviation of the average source value from the average of 
the values of several pixel pairs touching the pixel pair in 
question. Such test locates more accurately the boundaries of 
the region R. 

The most time-consuming step of the third method is the 
computation of the matrix V at all points. However, the ratio 
of its elements used in step (8) can be computed and 
tabulated beforehand as a two-argument function. In this case 
the computation of the matrix V is done quickly by 
interpolation and the whole noise reduction procedure takes 
about the same time as the calculation of image Ii from 
images Si, S 2 . 

The smooth fragments of the source images correspond to 
smooth fragments on reconstructed images (but not vice versa, 
since reconstructed images contain smaller number of details 
than the source ones) . Since the source images have 
relatively small noise, their use helps to define reliably 
the part of the region R, which does not contain details of 
the true image able to distort the average value, and helps 
to reduce statistical uncertainties of the boundaries. This 
makes the third method applicable in cases of very noisy 
(reconstructed) images, for which all other methods of noise 
reduction become inefficient. 

In cases of less noisy images, when different noise reduction 
methods are applicable, the third method makes smearing 
distortions in image I 3 smaller, and image I 4 cleaner from 
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shadows of fragment boundaries in image I 2 . These shadows may 
become even indistinguishable. 

The squared correlation coefficient for reconstructed images 
often exceeds 0.90 and noise dispersion can be reduced in 
image I 4 more than a dozen times, if matrix V is known 
accurately enough. The knowledge of reconstruction functions 
a, b gives the possibility to compute V and the ratio between 
fluctuations H and F more accurately than by analysis of 
images I lf I 2 themselves and achieve the noise reduction 
close to the theoretical limit. 

In an experiment of the third method, an image Ii and an 
image I 2 were .reconstructed from two source images of 
simulated X-ray pictures of the combination of bone and 
tissue. 

Reconstruction greatly amplifies the small intensity- 
dependent noise present in a source image (image Si or S 2 ) , 
so the noise in images Ii, I 2 was large and boundaries of 
fragments were smeared. This noise in Ii was 95% correlated 
with the noise in I 2 . In this case, difficult to the first 
method, the third method remained efficient and gave an image 
were the "shadows" of other image are noticeable, but the 
noise reduction is deep and small details are reproduced 
without any smearing. 

Preferably, the methods of the present invention is performed 
by software code, located in the internal memory of a 
computer, and executed by a processor of that computer. 

It will be obvious that the present invention may be varied 
in a plurality of ways. Such variations are not to be 
regarded as departure from the scope of the present 
invention. All such variations as would be obvious to one 
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skilled in the art are intended to be included within the 
scope of the appended claims. 



20 



CIAIMS 

1. A method for reduction of noise in an image including a 
plurality of pixels, comprising averaging pixel values over a 
region (R) , characterized by the steps of: 

5 - adding a selected pixel to the region (R) ; 

- grouping pixels adjacent the region (R) in pairs, wherein 
the pixels of each pair being oppositely located with respect 
to said selected pixel; 

- adding said pairs, pair by pair, to the region (R) in 
10 dependence on that the squared difference of the selected 

pixel value from ---the pairs half sums does not exceed the 
dispersion (D) of the noise of said difference multiplied by 
a tolerance level (L) ; 

- repeating said step of grouping and said step of adding 
15 said pairs until that, in said step of adding said pairs, the 

condition for adding said pairs is not fulfilled for any 
pair; 

- averaging the pixel values of said region (R) ; and 

- using the thus averaged pixel value for said selected pixel 
20 in reconstruction of said image. 

2. The method as claimed in claim 1, wherein said step of 
grouping excludes grouping of pixels previously being grouped 
in pairs. 

3. The method as claimed in claim 1, wherein said step of 
25 adding said pairs excludes pairs that do not touch any of the 

pairs already included in said region (R) . 

4. The method as claimed in claim 1, 2 or 3, wherein said 
method is performed a second time and wherein said step of 
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grouping during said second time only includes pixels 
rejected during the first performance of said method. 

5. The method as claimed in any of claims 1-4, wherein said 
step of adding said pairs is performed in dependence on that 
the squared difference of an average of pixel values in the 
region (R) from the pairs half sum does not exceed the 
dispersion (D) of said difference multiplied by a tolerance 
level (L) . 

6. A method for reduction of noise in an image including a 
plurality of pixels, characterized by the steps of: 

- obtaining a noise reduced value of a first pixel; 

- calculating an estimate of the fluctuation of said noise 
reduced value^of said first pixel; 

- calculating an estimate of the fluctuation at a second 
pixel, wherein said fluctuation of said second pixel is 
correlated to said fluctuation of said first pixel; and 

- obtaining a noise reduced value of said second pixel by 
subtracting said fluctuation at said second pixel. 

7. The method as claimed in claim 6, wherein the absolute 
value of said correlation is at least 0.8, preferably 0.9, 
and more preferably 0.95. 

8. The method as claimed in claim 6 or 7, wherein said first 
and said second pixels are located in different images. 

9. A method for reduction of noise in an image including a 
plurality of pixels, comprising averaging pixel values over a 
first region (R) around a selected pixel (X) , 
characterized by the steps of : 
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- finding a second pixel (XI) corresponding to said selected 
pixel (X) ; 

- adding said second pixel (XI) to a second region (Rl) ; 

- grouping pixels adjacent the second region (Rl) in pairs, 
5 wherein the pixels of each pair being oppositely located with 

respect to said selected pixel; 

- adding said pairs, pair by pair, to the second region (Rl) 
in dependence on that the squared difference of the selected 
pixel value from the pairs half sums does not exceed the 

10 dispersion (D) of the noise of said difference multiplied by 
a tolerance level (L) ; 

- repeating said step of grouping and said step of adding 
said pairs until that, in said step of adding said pairs, the 
condition for adding said pairs are not fulfilled for any 

15 pair; 

- averaging the pixel values of the first region (R) , which 
corresponds to the second region (Rl) ; and 

- using the thus averaged pixel value for the selected pixel 
(X) of said first region (R) in reconstruction of said image. 

20 10, The method as claimed in claim 9, wherein said step of 
grouping excludes pixels previously grouped in pairs. 

11. The method as claimed in claim 9 or 10, wherein said 
dispersion (D) is based on the region (Rl) instead of said 
selected pixel. 

25 12. The method as claimed in any of claims 9-11, wherein said 
step of adding said pairs is performed in dependence on that 
the squared difference of an average of pixel values in the 
region (Rl) from the pairs half sum does not exceed the 
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dispersion (D) of said difference multiplied by a tolerance 
level (L) . 

13. The method as claimed in any of claims 9-12, wherein said 
second region (Rl) is located in an image, which is different 

5 than the image wherein said first region (R) is located. 

14. The method as claimed in any of claims 9-13, wherein said 
noise reduced image is further noise reduced by the method 
according to any of claims 6-8. 

15. The method as claimed in any of claims 9-14, wherein said 
10 noise reduced image is further noise reduced by the method 

according to any of claims 1-5. 

16. A computer program product directly loadable into the 
internal memory of -a computer, said computer program product 
comprising software code portions for performing the method 

15 as claimed in any of claims 1-15 when said computer program 
product is run on said computer. 
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ABSTRACT 

The present invention relates to a method for reduction of 
noise in an image including a plurality of pixels, comprising 
averaging pixel values over a region (R) , comprising the 
5 steps of: adding a selected pixel to the region (R) ; grouping 
pixels adj acent the region (R) in pairs , wherein the pixels 
of each pair being oppositely located with respect to said 
selected pixel; adding said pairs, pair by pair, to the 
region (R) in dependence on the squared difference of the 

10 selected pixel value from the pairs half sums does not exceed 
the dispersion (D) of the noise of said difference multiplied 
by a tolerance level (L) ; repeating said step of grouping and 
said step of addi-ng said pairs until that, in said step of 
adding said pairs, the condition for adding said pairs is not 

15 fulfilled for any pair; averaging the pixel values of said 
region (R) ; and using the thus averaged pixel value for the 
selected pixel of said region (R) in reconstruction of said 
image . 
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